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Abstract 

Using Af = 2 superspace techniques we compute the four-loop spectrum 
^ ' of single trace operators in the SU{2) x SU{2) sector of ABJM and ABJ 

■ supersymmetric Chern-Simons theories. Our computation yields a four- 

loop contribution to the function /i^(A) (and its ABJ generalization) in the 
magnon dispersion relation which has fixed maximum transcendentality and 
coincides with the findings in components given in the revised versions of 
arXiv:0908.2463 and arXiv:0912.3460, We also discuss possible scenarios for 
an all-loop function h?{\) that interpolates between weak and strong cou- 
plings. 
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1 Introduction 



The ABJM model is an A/^ = 6 supersymmetric U{N) x U{N) Chern-Simons theory 
with opposite levels coupled to matter [1]. Like its cousin Af = 4 super Yang-Mills in 
four dimensions, its two point functions of single trace operators map to an integrable 
system in the planar limit [2H5]. For M = 4 SYM, the integrability has been used as a 
powerful tool to interpolate between strong and weak coupling, where one can see the 
perturbative behavior of the gauge theory morph into the stringy behavior expected 
from the AdS/CFT conjecture PE]. 

The ABJM model has two extra features that give it a richer structure than A/" = 4 
SYM, at least as far as the integrability of the two point functions is concerned. The 
first is that the Bethe equations and the dispersion relations contain an undetermined 
function /i^(A) of the 't Hooft coupling, A = N/k, where k is the Chern-Simons level [5]. 
The second is that the theory can be deformed into a U{M) x U{N) gauge theory while 
still maintaining the M = Q supersymmetry [8]. In this ABJ case there are now two 
't Hooft parameters, 

, M ^ AT , ^ 

and, if integrability is maintained, a single function /i^(A,cr), where 

A = \/aI, a=^^. (1.2) 

A 

The spin-chain that appears in the ABJ(M) models has 05*^(614) symmetry and 
is of alternating type, with the spins on the odd sites in the singleton representation 
of the supergroup and the spins on the even sites in the anti-singleton representation 
[2] 4ill9UT0] . In order to find h^{\, a) it is only necessary to consider the compact subgroup 
SU{2) X SU{2) of OS'p(6|4), with the spins on the odd sites transforming in the (2, 1) 
representation and the spins on the even sites transforming in the (1, 2) representation. 
The ground state has all spins aligned and the excitations (or magnons) are flipped 
spins that live on either odd or even sites. The dispersion relations for these two types 
of magnons are given by 



EoM = Vg2+4/l2(A,a)sin2|-g , Eeven(p) =^odd(p)L^_,, (1-3) 



where (5 = 1/2 for fundamental magnons while larger values of Q correspond to magnon 
bound states. 

At weak coupling the function /i^(A, cr) can be computed perturbatively. The leading 
contribution appears at two-loop order and is relatively easy to compute, both for ABJM 
dH, and ABJ dOlIII], where one finds 

hWa) = \^ + 0{y). (1.4) 

However, at strong coupling on the ABJM slice where a = 0, one readily finds from the 
string sigma model [3|[T2|[T3]. 

/.2(A,0) = ^A + O(l). (1.5) 
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Hence, h^{\,a) is an interpolating function and can be expected to have corrections at 
every even order of perturbation theory, with a general structure 



oo 



(1.6) 



n=2 




(1.7) 



This calculation was done using the explicit component action and involved the com- 
putation of dozens of Feynman diagrams. A straightforward extension of the methods 
in [14[IT5] to higher loops would lead to a mind boggling number of diagrams. Moreover, 
one would like to verify (or disprove) that the ABJ theory is integrable, even at the 
four-loop order. The SU{2) x SU{2) sector is trivially integrable at four loops, so it 
would be necessary to go beyond this sector to find a nontrivial check of integrability 
at this order. But even this seemingly modest task is extremely daunting in component 
language. 

In this paper we compute h^i^a) in (11. 7p using the superspace formalism. Superspace 
techniques have proven to be very effective in computing the dilatation operator [TB] 
and in evaluating wrapping corrections [TTIIIH] in A/" = 4 SYM [T9l|20] and in its (3- 
deformation [2TI423] . Naturally, one would also like to apply them to the ABJ(M) 
models. Their main virtue is that they drastically reduce the number of Feynman 
diagrams that one must compute. We will later summarize several restrictions on the 
allowed diagrams [16] that greatly limit the number that can contribute to h^{a). As 
we will see in this paper, at the two-loop order there is only one diagram in superspace 
that contributes to h^iX, cr). At the four-loop order there are 15 (plus reflections of some 
of the diagrams). Contrast this to the component calculation in [lll[l5], where one has 
many times more diagrams. Not only does this demonstrate the formalism's power, but 
it is also crucial in verifying that (11. 7p is actually correct (see footnote [T]) . 

One can also see from (I1.7P that hi{(y) has uniform transcendentality two. From 
the component point of view this seems almost miraculous since many diagrams have 
rational coefficients (that is, they have transcendentality zero), others have transcen- 
dentality two, and some are mixed. When everything is combined one finds that the 
rational coefficients cancel. In superspace, while there are still diagrams with rational 
coefficients, their cancellation appears more natural due to correlations between the 
single and double poles. 

We will also present two possible scenarios for an all-loop function for /i^(A), including 
one that might work. It reproduces the first two orders of perturbation theory as well 
as the leading sigma-model contribution at strong-coupling. The one-loop sigma-model 
contribution to /i^(A) depends on how a sum is carried out over an infinite number of 
modes. Our proposal disagrees with the more conventional prescription in [2l], but 
agrees with the prescription in [25]. The other proposal looks for a connection with 

different result for ft.4((T) in (|1.7p was given in earlier versions of [T31II1]. After it became clear 
that those results were in conflict with the results presented in this paper, an overall sign error was 
discovered for three of the Feynman graphs. 
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matrix models on a Lens space. These arise in the study of supersymmetric Wilson 
loops in ABJ(M) models [261429] . In particular, we consider the free energy of the 
matrix model which is a function of A. We will see that h'^{X) has a structure similar 
to the derivative of the matrix model free energy, both at small and large A. But the 
coefficients in their respective expansions do not quite line up. 

In order to complete the four-loop analysis in the SU{2) x SU{2) subsector, we 
will apply the superspace formalism to compute the leading wrapping corrections for a 
length four operator in the (1,1) representation of SU{2) x SU{2). Here we find that 
the wrapping corrections per se differ from those computed in component language. 
However, other range five interactions must be subtracted and this subtracted piece 
also differs from the corresponding term in the component calculation. The two effects 
combine to give the same four-loop anomalous dimension for this operator as was found 
using components. 

The rest of the paper is organized as follows: In section [2] we review the ABJ(M) 
models in A/" = 2 superspace. In section [3] we discuss the relation of the dilatation 
operator to h'^{X,a). In section H] we enumerate and compute all Feynman diagrams 
that contribute to the four- loop term h^^a). In section [5] we discuss our investigation 
into possible all-loop functions for h?{X). In section[6]we apply the superspace formalism 
to the wrapping corrections for operators of length four. In section [7] we present our 
conclusions, which includes suggestions for further work. Many further details, including 
the four-loop decoupling of odd and even site magnons and the consistency of double 
poles due to UV sub divergences can be found in the several appendices. 

2 ABJ(M) models in A/^ = 2 superspace 

In this section we review the Af = 2 superspace formulation for A/" = 6 superconformal 
Chern-Simons theory. This was first given in [30|, but in this paper we follow the 
notations used in [31j which are adapted from the ones of [32]. For the first papers on 
the M = 2 superspace formulation of Chern-Simons theory coupled to matter see [55H35] . 
Appendix [A] collects our notation and conventions. 

The U{M) X U{N) supersymmetric Chern-Simons theory has two M = 2 vector 
supermultiplets, V and V, with V transforming in the adjoint of U{M) and V in the 
adjoint of U{N). In order to extend the supersymmetry to A^ = 6, the ABJ(M) action 
also contains two sets of chiral matter superfields, and Wa with A = 1,2. and 
Wa transform respectively in the (2, 1) and (1, 2) of the global SU{2) x SU{2) fiavour 
subgroup described in the introduction. Moreover, they transform in the bifundamental 
representations (M, N) and (N, M) of the U (M) x U{N) gauge group. 

Each gauge group in the gauge fixed Af = 2 superspace action has associated with 
it a pair of chiral ghost superfields, c, c' for U{M) and c, c' for U{N) Including 
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all of these ingredients, the gauge fixed ABJ(M) action in A/" = 2 superspace reads 



5'cs + S, 



gf 



k 



J d^x d^e J^' dt tr V ( e~*^ D„ e*^ +^ +1 ) V 



d:'xd^e I dttrf fD^e-^^D^e^^+V^D^+iD^V 



5, 



FP 







>S'mat — 


-/ 

to 7 


'S'pot = 


k i ' 


4^2 . 



d^a; d^e tr(c' + c') L ly (c + c + coth Li^^(c - c)) 
d^x d^6' tr(c' + c') Li^(c + c + coth L i^(c - c)) 

2 2 

d^a; d^0 tr (Z^ e^ e"^ +1^^ e^ Wb ) 

d^a; d^^eAce^^ tr Z'^WbZ^Wd 
+ /" d^a;d¥e^^eBDtrZAiy^^cW^^ 



(2.1) 



where hy X = \V ,X] and a and d are gauge fixing parameters. 

Many of the terms in this action have an infinite expansion, but for our purposes it 
is only necessary to retain the first few orders of any expansion. The first term in the 
Chern-Simons Lagrangian expands to 



dttr1/D"e"*^D„e* 



itrVD^D^F-^trFD" [l^,D, V] + 
2 o 



(2.2) 



The quadratic piece in this expression, together with the a- and a-dependent gauge 
fixing terms, determines the gauge superfield propagators. In order to simplify the D- 
algebra manipulations we will choose the Landau gauge where a = d = 0. The leading 
expansion for the Fadeev-Popov action is 



S- 



FP 



k 

Air 



j d^x d'^e tr {c'c - dc + ^(c' + c') [l^,c + c] ) + 



(2.3) 



while the leading expansion for the matter action D-terms is 



5, 



mat 



k 

Att 



d^xdW 



tr Za{Z'^ + VZ^ - Z^V + -{V^Z^ + Z^V^) - VZ'^Vj + . 
+ trW''{WA + VWa - WaV + ]^{y^WA + WaV^) - VWaV 



+ . 



(2.4) 

We have collected the Feynman rules which follow from the action and the above 
expansions in appendix [Bl The supergraphs are then constructed from the Feynman 
rules and are reducible to ordinary integrals using standard D-algebra techniques [32] . 

The advantage of using superspace as opposed to the component approach is that the 
number of diagrams is significantly smaller. Furthermore, one can often find cancellation 
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patterns between different supergraphs or demonstrate finiteness theorems for classes 
of diagrams p^[T9] . Sucli generalized finiteness conditions [16] that follow from power 
counting arguments and some of their implications are summarized in section 14.11 They 
predict the finiteness of many diagrams and will be of great use to us in our calculations. 



3 The dilatation operator and /i^(A, a 



The dilatation operator T> is the natural tool to study the anomalous dimensions of 
composite operators in field theory. It can be defined as the operator that by acting on 
composite operators Oa provides the matrix of scaling dimensions 



(3.1) 



Note that A^'' leads in general to the mixing between operators. As known, the ma- 
trix of dimensions, and therefore the dilatation operator, can be extracted from the 
perturbative renormalization of the composite operators Oa 



Za'O,^ 



bare 



z = i + + x'Zi + 



(3.2) 



The matrix Z is such that Oa.ren is free from perturbative quantum divergences and can 
be computed in perturbation theory by means of standard methods. In this paper we 
use dimensional reduction with the space-time dimension D given by 



D = ^-2e 



(3.3) 



in order to regularize quantum divergences that show up as inverse powers of e in the 
limit £ — )■ 0. By introducing the 't Hooft mass and the dimensionful combination A/x^^ 
the dilatation operator is then extracted from Z as 



= I'ciassicai + yU^" In Z{Xfi e 
d/i 



^classical + hm 



2e\^\nZ{\e) 
dA 



(3.4) 



In a loop expansion of the dilatation operator, the /th loop order is then simply given 
by the A^' coefficient of the 1/e pole of hiZ multiplied by 21. The higher order poles 
must be absent in hiZ] this will be later used as a consistency check for our result. 

As discussed in the introduction, in the ABJ(M) models the dilatation operator can 
be mapped to the long range Hamiltonian of a spin-chain system for the whole 0Sp{Q\4) 
symmetry group [21IID]- We focus on the SU{2) x SU{2) subsector where the magnons 
propagating along the spin chain form two sectors: the ones living on the odd sites belong 
to the first SU{2), while those on the even sites are associated with the other SU{2). 
As demonstrated in appendix IF.lt in our four-loop analysis the two different types of 
magnons can be regarded as non-interacting, since the contributions to the dilatation 
operator of the respective diagrams that could lead to these interactions cancel. The 
all-loop Bethe Ansatz [5] predicts that such interactions start at eight loops. In analogy 
with the M' = 4 case, the spin-chain is interpreted as a quantum mechanical system in 
which the ground state of length 2L can be chosen to be 



(3.5) 
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With a single excitation W2 of an odd site the momentum eigenstate is defined as 

L-l 

= J2e'''\WiZ^fW2Z\WiZY-''-^ (3.6) 

This describes a single magnon excitation with momentum p. The main difference 
between the Af = 6 CS and the Af = 4 SYM case is the existence in the former of two 
different SU{2) excitations corresponding to the sectors mentioned above. 

Up to four loops, the dilatation operator for a chain of lenght 2L then expands as 

V = L + X\V2Md + V2,e.cn) + X^V^Mdi^r) + ^^4,even(fx)) + O(A^) , (3.7) 

where the individual parts act non-trivially on odd and even sites only. 

In the J\f = 4 SYM case chiral functions have been introduced in [19] as a very 
convenient basis for the dilatation operator of the SU{2) subsector. The chiral functions 
directly capture the structure of the chiral superfields in the Feynman diagrams. As in 
the A/" = 4 SYM case, also in the Af = 6 CS case the elementary building block for the 
chiral function of the SU{2) x SU{2) subsector is constructed from the superpotential 
by contracting one chiral and one anti-chiral vertex with a single chiral propagator. The 
resulting flavour structure then yields the simplest non-trivial chiral function. 

The chiral functions that are relevant to two loops in A/" = 4 SYM and to four 
loops in A/" = 6 CS theory turn out to have identical form in terms of the respective 
permutation structures and read 

X(a,6) = {a,b} - {a} - {b} + {} , 

X(«) = {a} - {} , (3.8) 

XO = {} . 

However, the permutation structures in both theories slightly differ. In the Af = 6 CS 
case they are given by [15] 

L-l 

{ai, a2, . . . , dm} = P2«+ai 2i+ai+2 ^21+02 2i+a2+2 ■ ■ ■ P2i+a,n 2i+a™,+2 ; (3-9) 

i=0 

where we identify L + i = i, such that the product of permutations, in which Paa+2 
permutes the flavours at sites a and a + 2, is inserted at every second site of the cychc 
spin chain of length 2L^ The insertion at each second site thereby allows for the 
decomposition of the dilatation operator into two separate pieces acting only on odd or 
even sites as in f l3.7p . The decomposition of the dilatation operator to four loops [TJ] in 
terms of chiral functions then reads 

'C2,odd = -x(l) , 

^^4,odd(a) = -x(l, 3) - x(3, 1) + (2 - K{a))x{l) , ^ ' ' 

I^4,even(cr) = -x(2, 4) - x(4, 2) + (2 - h,i-a))xi2) . 



^Note that the permutation structures obey 

{...,a,6,...} = {...,&, a,... } , \a-b\^2 



{a, . . . ,b} ^ {a + 2n, . . . ,b + 2n} . ^^'^^^ 
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The coefficients are thereby fixed by the magnon dispersion relation (11 .Sp in terms of 
the four-loop contribution h^i^cr) of the a priori undetermined function /i^(A, a) in (11. 6p . 
As explained in [15] to obtain the above result, one just has to compare the expansion 
of the magnon dispersion relation to the momentum dependence when the individual 
terms are applied to the single magnon momentum eigenstate (I3.6p . 

The function h4{a) can be computed in the weak coupling limit from a direct per- 
turbative calculation. This has been done by using component fields techniques in [T5] . 
Here we present its calculation by using M = 2 supergraphs. As in the component 
calculation [15], also here it suffices to only consider the odd part of the dilatation oper- 
ator, i.e. the contributions with chiral functions that have odd integers as argumentsE] 
The supergraphs computation of the full r'4,odd, and in particular of h^i^a), is the main 
result of our paper. 

4 Feynman diagram calculation 

Before starting with the explicit evaluation of Feynman diagrams we will summarize the 
previously mentioned finiteness conditions which allow us to disregard entire classes of 
diagrams. 

4.1 Finiteness conditions 

Based on power counting and structural properties of the Feynman rules, in [16] finite- 
ness conditions for Feynman diagrams of = 4 SYM theory in terms of A/" = 1 
superfields and for A/" = 6 CS theory in terms of M = 2 superfields were derived. They 
hold for each diagram that contributes to the renormalization of chiral operators in the 
respective SU (2) or SU (2) x SU (2) subsectors. In Landau gauge, such a diagram with 
interaction range R > 2 has no overall UV divergence, if at least one of the following 
criteria is matched^ 

1. All of its chiral vertices are part of any loop. 

2. One of its spinor derivative is brought outside the loops. 

3. The number of its spinor derivatives brought outside loops becomes equal or 
bigger than twice the number of chiral vertices that are not part of any loop. 

In the fiavour SU{2) x SU{2) subsector, a chiral vertex that is not part of any loop 
always generates fiavour permutations and therefore a non-trivial chiral structure of the 
diagram. Analogously to the A/" = 4 SYM case, the above finiteness conditions hence 
imply the following rule: 

• All diagrams with interaction range R >2 and trivial chiral structure xO are finite. 

^As we mentioned before, odd and even site magnons are decoupled here, there is therefore no 
contribution with chiral functions with both odd and even integer arguments. We explicitly demonstrate 
their absence at four loops in appendix lF.il 

^i? > 2 means, the composite operator is IPI connected with the rest of the diagram, not including 
the non-interacting fields of the operator. 
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Together with the conformal invariance on the quantum level, i.e. the finiteness of the 
chiral self energy, this implies that any diagram which does not manipulate the flavour, 
i.e. it has trivial chiral structure xO defined in (13.81) . has no overall UV divergence. 

Since the propagators of the vector fields in Landau gauge carry D D, the finiteness 
conditions imply the following statement: 

• A diagram with interaction range R>2 has no overall UV divergence, if it contains 
at least one cubic gauge-matter interaction with a chiral field line which is not part 
of any loop. In particular, if in the diagram exactly one of the chiral vertices appears 
outside the loops, then it also has no overall UV divergence if the anti-chiral field 
of at least one cubic gauge-matter interaction is not part of any loop. 

According to this statement, there are no contributions to the dilatation operator that 
come from diagrams in which the chiral line of a cubic gauge-matter vertex is an external 
line. In appendix [E] we will, however, evaluate such diagrams with IR divergences explic- 
itly to show that indeed all IR divergences cancel out in the renormalization constant 
Z in 



4.2 Two loops 

Before attacking the more involved four-loop case, let us see how the two-loop result is 
obtained from supergraphs. There is only one non-vanishing logarithmically divergent 
diagram contributing. It evaluates to 



{4:71 



XXI 



k 



2 MNhx{l) = -r-xa] 
^ A e 



(4.1) 



where the two-loop integral I2 is given in flC.2|) . As already discussed, to obtain the 
contribution to the dilatation operator one has to take the coefficient of the pole 1/e 
and multiply it by —21, in this case equal to —4. Once a factor = A A is removed one 

gets 

V2 = • (4.2) 



This coincides with the results found in 



in components. 



4.3 Four loops 

Now, let us move to the four-loop contributions to the dilatation operator. We will 
separate them according to the range of the interactions. We will explicitly present 
only the diagrams surviving the finiteness conditions of [16] that are summarized in 
section 14.11 It is important to note that, according to these arguments, an overall UV 
divergence can be present in superficially logarithmically divergent diagrams if at least 
one purely chiral vertex remains outside the loops. This implies that the minimum range 
of interaction at any loop is three. This is consistent with the fact that the minimal 
structure that can appear in the dilatation operator is The range varies between 

three and the maximum one which at four loops is five. 
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Note that together with the 1/e poles we will also keep the higher order poles that 
display the presence of sub divergences. Here, to four-loop order the only appearing 
higher order poles are double poles. In appendix IF. 21 their cancellation in In Z will be 
explicitly demonstrated as an important consistency check of our calculation. 

We note that, for the convenience of the reader, all the integrals appearing in the 
following are collected in the appendix O 



4.3.1 Range five interactions 

At four loops there is only one supergraph that involves the maximum number of five 
neighbouring fields in the interaction. It is given by 

- ^(MiV)'/a(1.3) =(§!(- ^ + H),(l,3) . (4.3) 

By taking into account the reflected diagram, the maximum range contribution to the 
renormalization constant i^ 

2r5,odd = -(1 + 7^)Sr5 = ^ - ^) 3) + x(3, 1)) . (4-4) 




4.3.2 Range four interactions 

There are four diagrams which have range four interactions and contribute to the struc- 
ture in the dilatation operator. According to section for an overall UV diver- 
gence to be present, at least one purely chiral vertex has to remain outside the loops, 
and a single gauge propagator can not end up on an external leg. Therefore, the only 
relevant contributions turn out to be 



'r4 



Vr 



r41 



Vr 



r42 



Vr 



r43 




M=^Ar/4bbbX(i: 



16 



Is 



32 



32 



(4.5) 



A^A / 1 1 7r2 



16 l2^-il2-T))^W 



^By 7?, we indicate the reflection of a supergraph at the vertical axis. As in [T^, the operation 
preserves the type of chiral function, i.e. if it belongs to the odd or even sector. In case of an even number 
of neighbours interacting with each other the operation therefore involves a shift of the interaction by 
one site along the composite operator. Effectively, Tl, therefore exchanges A with A and x(fl: ^) with 
a)- 
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Also in this case one has to consider the diagrams obtained by reflecting the previous 
ones. The total contribution to the renormalization constant is then 



4.3.3 Range three interactions 

The range three interactions arise from two-loop corrections to the propagators and 
vertices involved in the two-loop diagram f l4.ip . It is important to note that, due to the 
finiteness rules of section I4.H overall UV divergences can arise only from corrections to 
the lower vertex or one of the three lower chiral propagators. According to the analysis 
of [31], the two- loop corrections to the chiral two- and four-point functions are plagued 
by IR divergences even if free of UV poles. This is due to the particular structure of 
the gauge superfield propagator and cubic vertices in A/" = 2 superspace. We stress 
that IR divergences do not appear in component fields [15], since in three dimensions 
IR dangerous cubic vertices contribute non-trivial momentum factors to the numerators 
of the loop integrals. In superspace, the appearance of IR divergences in intermediate 
steps can be cured by using a non-standard gauge fixing procedure first introduced 
in four dimensions in [10] and adapted in [31] to the three-dimensional case. Since 
we are interested only in the overall UV divergences of the diagrams, a computational 
strategy could be to ignore purely IR divergent diagrams and to IR-regulate diagrams 
that involve both UV and IR divergences in such a way as to extract the purely UV 
poles. For example, this is illustrated in appendix IC.2I where we can regulate the IR 
divergences by inserting external momenta in IR divergent diagrams. However, in the 
main body of the paper we have decided to keep track of the IR divergences and check 
at the end their cancellation. Such a check is described in appendix |El 

The interested reader should look at appendix [D] for a description of the two-loop 
corrections to the two- and four-point functions needed in the calculations of this section. 




(4.6) 
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The contributions with only UV divergences are given by 

My 



'r3 



r316 



r32a 



Vr 



r326 



r34 



K^31a — 




2k^ 



1 1\ . ^ 



/4\ 




K-33a — 



(4 



16 



M3iV(/4 + /42bM)x(l) 



16 4£ 



x(i), 



(4 



16 



(4.7) 



(MiV)2 J422qtrAS6dX(l) 



1 1/^ 27r2 



))x(i) , 



^336 — 



^ j^} (-^-^)^(2/42bbe — -^422qtrAB6(i 



(AA)2 / TT- 



16 V 3£ 



+ 2(2/22i6e - /22idc)G'(2 - 2A, 1)G(2 - 3A, 1) 

— 2(/42bbd + -^42b6e))x(l) 

x(i). 
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The contributions with both UV and IR divergences are given by 



Vr 



r35 



Vr 



My 



{MN{4MN - M^)) [h - /4UVIR + /42b6d)x(l) 



AA 



r36 



16 



(4AA - A^) 



1 

2^2 



TT 



7-ln47r W(l) 



MiV(^2MiV/4bbb - -(8MiV - (M^ + iV2))/4uviRjx(l) 



= (^^T + ('^^ - (A^ + A^)) + ^ (2 - 7 + In 4.) ) ) x(l) . 

(4.8) 

Note that the expressions for the integrals that appear in the results have their UV 
sub divergences subtracted. The suffix UVIR appears on integrals which due to different 
arrangements of their external momenta contribute both UV and IR divergences. The 
UV poles can be extracted by adding external momentum to the cubic vertex which 
causes the IR divergence, i.e. one replaces /4UVIR — ^ -^4- This then yields 



UV 

r35 



^(4AA-A^)f--L + if2-^))x(l) 



UV 
r36 



16 

AA / - TT^ 

— AA h 

16 V e 



2e' 



4 



m-(A^ + P))(^-i 



(4.9) 



In appendix |E] we explicitly demonstrate that this result is also obtained if instead of 
choosing an IR safe momentum configuration all relevant diagrams with IR divergence 
are considered, i.e. the IR divergences cancel out in the final result. 

The contribution of the range three interactions to the renormalization constant Z 
is then given by 



rS.odd 



-(1+ 7^)(5r3 + K31a + l^r316 + V^a + 2^326 + 2^34 + K^^) 
— K33a — Vr33b 



OT/UV 



AA 

16 



AA 



1 1 / 27r2 
- + - 4+ — 

e\ 3 



(4.10) 



TT 



H-(A^ + A%)X(1) 



4.4 Final result 



We are now ready to put together the parts of our calculations necessary to extract the 
four-loop dilatation operator. As discussed before the dilatation operator for odd sites 
is obtained by extracting the 1/e pole from the renormalization constant. Summing up 
the contributions to the 1/e pole from (14.41) . (14.61) and (14.101) . we obtain 



A^-Z4,odd I i 



(^r5,odd "I" ^r4,odd 

AA 



2, 



r3,odd 



)l 



165 



2AA(x(l,3) + x(3,l))+ AA 4 + 



27r^ 



+ (A^ + A^)- X(l 



(4.11) 
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that, rewritten in terms of A and a of (11.21) . gives 

. , , ^2. ^2. 



A -Z4 odd 1 ~ T7^ 



- 2(x(l, 3) + x(3, 1)) + (4(1 + ^) + ^'y • (4-12) 



As already observed, in the ln.Z the higher order poles must be absent. This is a 
useful consistency check of our computation. Additional diagrams that do not contribute 
to the dilatation operator but have non-vanishing double poles have to be taken into 
account. Some of them consist of two separate two-loop interactions. Furthermore, one 
has to consider the diagrams that lead to interactions between magnons at odd and even 
sites and contribute only to the double pole when summed up. In appendix [0 we prove 
that when all these double poles are taken into account, their sum is indeed cancelled 
by the two-loop contribution in the expansion of In^. The dilatation operator for odd 
sites is then obtained from (14.121) by multiplying the 1/e pole by 8. With ^(2) = it 
reads 

I^4,odd(fx) = (2 + (4 + a2)C(2))x(l) - 3) - x(3, 1) . (4.13) 

By comparing the previous result with equation (13. lip we read off the four-loop coeffi- 
cient of the function /i^(A, a) 

/i4(a) = -(4 + a^)C(2) . (4.14) 

This result coincides with the one computed in [15]. It is interesting to note that, in 
contrast to the component calculation in [15], the integrals that contribute here to the 
dilatation operator show a correlation between the quadratic and the rational simple 
pole in e: their relative coefficient is always —4 as for the simplest four-loop integral 
1^ in (lC.4p . The rational term in (I4.13P and therefore its absence in (I4.14p is hence 



correlated with the quadratic pole that itself is determined by the two- loop result ( ]F.5p . 



5 Possible scenarios for an all-loop function 

In this section we discuss our attempts to find an all-loop function for /i^(A, a). 

In the ABJM case where cr = 0, /i^(A,0) = h'^{\), there is a surprisingly simple 
function that matches the weak coupling behavior up to four-loop order and also matches 
the leading strong coupling behavior. To this end we define t = 27rzA, which is a 
natural variable that also appears in expressions for supersymmetric ABJ(M) Wilson 
loops [26l[28l|29] . We then consider a rescaled function g{t) = {271)"^ h? {\) . In terms of 
g{t) the magnon dispersion relation becomes 

and so has a form more in line with the TV = 4 dispersion relation where in that case 
g{t) in (15. ip is replaced with A. 

In terms of g{t), the proposed all- loop function is 

g{t) = -(1 - t) log(l -t)-{l + t) log(l + t) , (5.2) 
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whose weak coupling expansion is 



l^^n{2n-l) 6 15 ^ ' (5.3) 

= (27r)2 (A' - 4 C(2) + 6 C(4) A' + C(A')) . 

An obvious test is to compute /i^(A) to six-loop order, where the all-loop function in 
(15. 2 p predicts the value = -^^f^. A six-loop computation is admittedly very difficult, 
but we believe it is manageable using the Af = 2 superspace formulation. 
At strong coupling the expansion is 

g{t) = -ivrt - 2 logt - 2 + 0{r^) 

(5.4) 



(2^)' - (^ l°g^27rA) - 2 + 0{X-')^ . 



The dominant term agrees with the leading strong coupling expansion from the string 
sigma-model. But also observe that the ffist correction corresponds to a two-loop con- 
tribution; a one-loop correction is absent. This disagrees with the prediction in [2l] 
arising from the one-loop correction to the energy for a folded-string [25 |I^TI - H^ . In this 
language one would expect a g(t) with leading asymptotic expansion 



g(t) = -mt- 2 V^iTrt ln(2) + . . . . (5.5) 

However, if one chooses a different prescription for summing over mode frequencies, 
where one essentially groups the modes into heavy and light [25j, then g{t) no longer 
has the -y/t term, agreeing with the large t expansion f l5.4p rl 

The function in (15. 2p does not appear to have an easy generalization to the ABJ case 
where cr 7^ 0. Such a function would be expected to be invariant under the transforma- 
tion [g 

A^A, A^2A-A + 1. (5.6) 

Under (15. 6 p the perturbative regime is mapped into strong coupling, making its ver- 
ification difficult. Some evidence that /i^(A,cr) is consistent with (15. 6p was presented 
in IIH]. One possible hint about the all- loop structure is that the four-loop contribution 
to h'^(X,a) can be rewritten as 

+ a^) = XX{X + (5.7) 

which is zero if A = —A It would be interesting to see if the higher order corrections 
remain zero under this condition. However, it is not clear how this could square with 
the strong coupling behavior nor with an invariance under the transformation in (15. 6p . 

Another possibility is that h?{X,(j) is somehow related to recent results concerning 
supersymmetric Wilson loops in the ABJ(M) models. In this latter case, it was found 
using localization [191150] that the Wilson loop expectation value could be reduced to 



^See [15] for a further discussion of this. These authors also show that the same choices of pre- 
scriptions appear in finite size corrections for giant magnons |46[|47| and lead to the same one-loop 
contributions to h^{X). 
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a matrix model on a Lens space [26]. This matrix model is solvable in the planar 
limit [5T1[52] and hence all-loop predictions can be extracted. In particular, for ABJM 
the perturbative free energy of the matrix model is |28] 

F{t) = N' (^log(t) + + Oit')^ . (5.8) 

It is tempting to look for a connection between F{t) and g{t). One might try 

The full expansion also is maximally transcendental, but here one finds that the term 
is off by a factor of 2/3. At strong coupling the free energy is asymptotically [29] 

F{t) ^ -AT^— - {-it)-^'^ . (5.10) 
Applying the same rule as in (15. 9p one finds 

(^w)^/^ = - f (—^)^/^ (5.11) 

which differs by an overall factor of 7r/3 from the square root of the leading term in 

dsai). 



6 Wrapping interactions 

To obtain the complete four-loop spectrum of operators in the SU (2) x SU (2) subsector, 
we have to consider the wrapping interactions for the non-protected operators that 
consist of up to four elementary fields. The only non-trivial operator is in the 20 of 
SU{A) and has L = 2, i.e. exactly four elementary fields. 

The only wrapping diagrams which according to the initially discussed finiteness 
theorems based on power counting can contribute to the dilatation operator are given 
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by 



Wo 



Wa 




^-^(MiVfJ4x(l) 



2(47r) 



16 Ve^ £ 



^ ^^(MiV)2j422btrABcdX(l) 



16 



16 



2(47r) 



-(MiV)2/4x(l) 



422qtrA_BbdXUj — — TTT" 
16 



If 2 , 

-(4- -TT^ 

£ V 3 



x(i) 



(6.1) 

There are four distinct diagrams of type W2 and two of type W^. The sum of the 
wrapping diagrams is therefore given by 
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7r 



xii)- 



(6.2) 



Multiplying the 1/e pole of by —8, we obtain the wrapping contribution to the 
dilatation operator. It reads 



V 



4,odd 



-(2-2C(2))x(l). 



(6.3) 



Now, by subtracting from fl4.13p the range five contribution and inserting hi{(T) 
— (4 + (7^)^(2), the subtracted dilatation operator becomes 



V 



sub 
4,odd 



/.4(a))x(l) = (2 + (4 + a2)C(2))x(l) 



The dilatation operator for length four states then reads 



V 



range 4/ 
4,odd \ 



^) = ^^rdd(^) + ^^4"odd = (6 + ^^)C(2)x(l) 



(6.4) 



(6.5) 



and it coincides with the results obtained in terms of component fields p^[T5] . 

Note that the separation of the dilatation operator into wrapping and subtracted 
parts differs in the superfield calculation from the one obtained in component fields 
in [TUfTS] . The sum of the two terms is, however, the same in the two calculations, and 
hence the resulting anomalous dimensions for operators with length 2L = 4 agree. 
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7 Conclusions 



In this paper we have computed h^la) using the M = 2 superspace formahsm. The com- 
putation is greatly simphfied from the component version [131 [15] because the manifest 
supersymmetry in combination with finiteness conditions leads to a large reduction in 
the number of Feynman diagrams. 

With this reduction in diagrams, it should be possible to tackle more challenging 
computations, including the six-loop term hQ^a). Six loops would give one more data 
point and might provide further insights into an all-loop function. 

Alternatively, one could also apply the superspace formalism to four loops but beyond 
the SU{2) X SU{2) sector. This would not give us further information on /i^(A, a), but it 
would provide a check of higher-loop integrability in both AB JM and AB J models. One 
reason that integrability in the ABJ case is not assured is because at strong coupling a 
nonzero a would correspond to a nonzero ^- angle for the world-sheet, which is normally 
thought to destroy integrability. However, at the lowest order in perturbation theory, the 
spin-chain is integrable in all sectors, even when cr 7^ [10l[TT]. It would be interesting 
to see how this plays out at higher loops. 
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A Conventions and identities 

We use three-dimensional spinor and superspace notations adapted from [32]. We 
directly work in the Wick rotated Euclidean space-time with metric g^y = g^^ = 
diag(l,l,l). For a given three-dimensional spinor field ipa, we raise and lower spinor 
indices as 

= C'l^^^ , ^„ = ^^Cp^ . (A.l) 
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where we use the spinor metric C"^ defined by 



For the contraction of spinor indices we use the notation 

= = X^V^a = XV^ , = \r'^a . (A.3) 

The 7-matrices obey the relation 

{iT.iiyp = - ^'""{i.Tp ■ (A.4) 

where the Levi-Civita tensor is such that e°^^ = 1. When one spinor index is lowered or 
raised the 7-matrices are symmetric 

= {Y)JCsp = , hT^ = c"\Y)s^ = {Yf"- . (A.5) 

The trace of product of 7-matrices satisfies 

tr(7V) = hT^iYfa = -23^" , 
tr(7'^7V) = -{iTpW^iYVa = -26^"' , (A.6) 
tr(7^V7'^) = {iTp{rY,{l''V5{l''Yc. = 2((7'^V^ - g'-'g'^ + g^'^gn . 

We use the convention that the first of two contracted indices is always an upper index; 
this is used in the previous formulas in the definition of the trace of products of gamma 
matrices and it is very useful for D-algebra manipulations [32] . 

Using the 7-matrices we can move from vector to bi-spinor indices thanks to the 
following definitions 



Pap = {Y)apP^. , = ^ilt^T^PalB , (A.7) 

respectively for coordinates, momenta and fields. As usual, here the momentum is 
related to the vector derivative 9^ = ^ by Fourier transform and = id^. 

The three-dimensional, J\f = 2 superspace spinor covariant derivatives Dq,, Dq, satisfy 
the algebra 

{D,,D4 = {D„,D4 = 0, {D,,D4=p„^. (A.8) 
The metric e^B for the SU(2) fiavour indices is given by 

ei2 = 1 , e'' = l, e^^ecD = 6^5^ - 5^5^ • (A.9) 
The fiavour indices are raised and lowered as 

^A^^AB^^^ ^A=^%A. (A. 10) 
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For the integration over the superspace our conventions are J — ^d^da, J — 
Id^'da emdjd^e = Jd^Od'^e, such that 

[ d^xd'^9 = [ d^x , f d^xd'^e = ! d^xIi\^e^Q , 

r J J (A.ll) 

d^xd^e^ / d^x D^&le^g^o ■ 



The 6'-space 5-function is given by 

5\e - e') ^{e- e'f{e - e'f . (A.12) 



B Feynman rules in superspace 

Wc use the Wick rotated Feynman rules, i.e. we have e"*"^ in the path integraL 

The propagators are given by 

{V{p)V{-p)) = -{V{p)V{-p)) = \BBS\e, - 62) , 



— - = {z^{p)z4-p)) = {w^{p)Wa{-p)) = ^^5\e, - 9,) , 



= (c-'(p)c(-p)) = -(c'(p)c-(-p)) 
= - {c'{p)c{-p)) = {c'{p)t{-p)) = ^^5\9, - 62) , 

where diagonahty in the gauge group indices and a factor ^ for each propagator have 
been suppressed. 

The vertices are obtained by taking the functional derivatives of the Wick rotated 
action (no factors of i) w.r.t. the corresponding superfields; we will give only the ver- 
tices involved in the computations of our paper. When a functional derivatives w.r.t. 
the (anti)-chiral superfields is taken, factors of (D^) are generated in the vertices. 
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Omitting factors for the three point vertices we obtain 





tr (T" [T^T'=] ) , 



yyzBzc 





VZbZC 




Vvcc' — 



Vvcc' — 



5: i tr (T-^ [T^T^] ) , Ws5' = 



5^ii{T'B,B^) , 

{-l)6^tr{T^B'-B^) , 




i tr (r« ) , 

5'' 1 tr (T'^ [T\T^] ) , 



(B.2) 

where the colour indices are labeled (a, b, c) counter clockwise starting with the leg to 
the left. Besides the matrices and transforming in the adjoint of the respective 
gauge groups U (M) and U {N), we have introduced matrices B- and Ba, with underlined 
a — 1, • • • , MN indices that transform in the (M,N) and (N, M) of the gauge group 
U{M) X U{N). The previous notations arc useful because one can effectively consider 
all the matrices to be the same for M = N and then only at the end one can easily 
recover the different factors of M and N coming from the colour contractions. 
The quartic vertices used in the paper are 



V^ZCZu — 




V^ZcZO 



^ U^[ir {{T\T']B^-Ba)] , 




(B.3) 



VZBVZ 



D <) 




-l)5^tr (T"S^™d) , 



where the colour indices are labeled (a, 6, c, d) counter clockwise starting with the leg 
in the upper left corner. The vertices Vj>2^y^^yD, Vyiwcwu'' ^vweVW^ involving the Wa 
and W'^ superfields are respectively identical to the previous three vertices up to trivial 
modifications in the flavour and colour structures. 
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The quartic superpotential vertices are 



^Z^WbZCWi 



D 




tr {B'-B,B^B^] , 



(B.4) 




teAce''''[ti {BaB'-B,B^) 



tr {B,B'-BaB^)] , 



where again the colour indices are labeled (a, b, c, d) counter clockwise starting with the 
leg in the upper left corner. Note also that, in a standard way, one of the (D^) factors 
has been absorbed into the (anti)chiral integration such that the integration measure of 
the (anti)chiral vertex is promoted to the full superspace measure. 

C Integrals 

In this section we collect the integrals required for our paper. The results are based 
on the Appendices H, I, J of [15] where the reader should look to have a complete 
description of the notations and results that we are using. 

The integrals are computed by using dimensional regularization in Euclidean space 
with D dimensions and 



As usual we will expand the integrals in the limit e — )■ up to the order needed for 
our computations. The parameter A in this appendix should not be confused with 
the 't Hooft coupling that appears in the main body of the paper. The integrals have a 
simple dependence on the external momentum which we will omit. Relations between 
four-loop expressions are understood to hold for the pole parts up to disregarded finite 
contributions. 

C.l Integrals with only UV divergences 

We need the following two-loop integral 



The reader can look at the appendix H of [15] for our notations in using the G-functions. 
Furthermore, we need the following two-loop integrals with two contracted momenta in 



D = 2(A + 1) = 3 - 2e , 




(C.l) 




(C.2) 
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their numerators 

W= =^(-Gi(l,l)G(l,l)-G(l,l)Gi(2-A,l) + Gi(l,l)Gi(2-A,l)) , 
W= =-Gi(l,l)Gi(2-A,l) . 

(C.3) 

At four loops there are many integrals involved in the computations. Here we list 
the results for the pole parts of the UV logarithmically divergent integrals where the 
subdivergences have already been subtracted. Four-loop integrals with no momenta in 
their numerators are 




Four-loop integrals with two contracted momenta in their numerators are 




Let us consider now four-loop integrals with four pairwise contracted momenta in their 
numerators. The following ones 



h22hABcd — 




(C.6) 
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appear in a fixed combination which can be recast into the form 



U22btrABcd 




422h ABcd ~ 422h AcBd + I i22h AdBc 




(C.7) 



(87r)4 V£2 



Here we have taken the trace of 7-matrices contracted with the momenta in the integraL 
We thereby read off the momenta in a cycle, but keep their direction as indicated by 
the arrows. 

We also need the integrals 



'422qABbd 



'422qAdB6 



lA22ciAbBd 




1 Y.i 



1/1 1 /5 vr 



iTry V4£2 £ V4 12 



(C.8) 



The linear combinations of integrals originating from the traces of 7-matrices read 



'422qtrAB6d 



= -tr (-e 



h22qtiABCD — tr 




-2(/422qAB6rf ~ h22ciAbBd + -^422qAdB6) 



— -(a -I 2 



1 TT^ 

(87r)4 3£ ■ 



There is an interesting relation involving the traces. It reads 



lA22qtxABCD — -^422qtrAB6d + 2/4 + ^lA2hbd 



1 Tl' 

^7r)4 2>e 



(C.9) 



(C.IO) 
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C.2 Integrals with IR divergences 

In this subsection we collect the integrals having poles in e which are due to IR diver- 
gences. By suffixes IR and UVIR we thereby label integrals which have one or both IR 
and UV divergences. 

The simplest two-loop integral with both an IR and an UV divergence is the loga- 
rithmically divergent tadpole 



'2tp — -'2UVIR 



. (C.ll) 



It is zero in dimensional regularization, i.e. the IR and the UV divergence cancel against 
each other. The UV divergence can be extracted by reshuffiing the external momentum. 
In particular, the UV divergence of /2tp is I2 defined in flC.2p : then the IR divergence of 
/2tp is -h- 

The simplest two-loop integral with only an IR divergence is given by 

1 / 1 



'2IR 




G(l, 1)G(2 - A, 1) = — ^ ( _ - + 2(1 + 7 - ln47r) + 0{e) 



One four-loop integral with both, an IR and a UV divergence is given bjlll 

WiR= fj\ =K{G{l,lfG{l-\l)G{l-2\,2-\)) 

^(-i + ^(-2 + 7-ln4.)). 



(C.12) 




(C.13) 



(87r)n 2e2 5 
Its IR divergence is extracted as 



'4UVIR 



where 1^ removes the overall UV divergence, since /4UVIR does not have a UV subdiver- 
gence. 

The simplest four-loop integral with only an IR divergence as overall divergence is 
given by 




■■■■R=-K )- =K(G(l,l)-G(l-A,l)G(2-2A,2-A))-K(/,)/2iH (^15) 

Here we have subtracted the UV sub divergence. 

^Note that, according to [15], with K() we mean the extraction of the pole parts of a function of e. 
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D Relevant one- and two- loop sub diagrams 

In this appendix we collect the results for the planar contributions to the one-loop 
vector superfield two-point function and the chiral superfield two-loop contributions to 
the two and four-point functions. The two-point functions have been first computed 
in [39l[53l[5l] in the Landau gauge. Such analysis has been extended to general gauges 
in [31] where the four point functions have also first been given for the ABJM case. Here 
we give the results extended to the U (M) x U (N) AB J case in the Landau gauge. The 
two-loop corrections to the chiral propagator and superpotential enter as subdiagrams 
in the evaluation of the dilatation operator given in section [31 

D.l One- loop vector two-point function 

For the U{M) vector superfield V the one-loop two-point function gets contributions 
from three kind of diagrams respectively having matter, ghosts and vector superfields 
propagating in the one-loop bubble. 

The contribution coming from the chiral matter superfields is 



The diagrams involving a loop of vectors sum up to the following contribution 




(D.l) 



The ghosts correction is 



^y,ghosti 




y, vectors 




The total contribution to the two-point function for the V superfield is then 




^ -5'^^G(1, 1)(4A^ - M) D" D, 



(D.4) 



The corrections to the U{N) gauge vector V two point function are clearly the same 
with the only difference that one has to exchange M with in the results. 
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D.2 Two-loop chiral two-point function 

The non-vanishing contributions to the two-point function of chiral superfields can be 
seen to arise from the following diagrams 

^ 2MNI2 , 
^ 2MNI2 , 

(D.5) 

^ -MiV(G(l,l))2 , 
^ (4iV - M)MG{1, l)Gi{l, 2 - A) 

= ^(4iV-M)M(/2tp-/2 + /2IR) 

^-^(4iV-M)MJ2tp , 

where, in each contribution, we have omitted a factor together with the colour 

and flavour structures. As discussed in section IC.2t the tadpole integral /2tp is zero in 
dimensional regularization. However, we keep track of it by splitting its UV and IR 
divergent parts. This is necessary for the check of the cancellation of the IR divergences 
performed in appendix |El 

Taking into account reflections of the diagrams at the vertical and horizontal axes 
where necessary, and summing up the contributions, the result reads 

Ec= • ^-2MN{G{l,l)f + h8MN-{M^ + N^))l2m. (D.6) 

Note that the result is UV finite and it includes an IR divergent term which turns out 
to be gauge dependent [3l] and, according to the discussion in appendix [El does not 
contribute to the dilatation operator. 
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D.3 Two-loop chiral four-point function 



The two-loop renormalization of the superpotential has been studied in [3T]. Here we 
summarize the results and extend them to the ABJ U{M) x U{N) case. It holds 



^-(47r)^A^(pi+P2) 




I 



(47r)2A2 




^ (47r)2AAtr(7'^7^7V) 




— > 



(47r)2AA tr(7'^7 V7 Vt^) 




tr(7'^7'^7"7^) 



, (D.7) 



Here the external momenta (pi, ■ ■ ■ ,^4) are ordered counterclockwise with pi the mo- 
mentum of the upper-left leg. 

The last contribution is rather complicated. However, it can be simplified by using 
momentum conservation to eliminate P2 the trace and the symmetrization inside the 
trace as 

^(tr(7^^7V7^^) + tr(7^^7^"7^")) ^ 
= -c/^'^tr(7V7V) + ^^''tr(7V7V) tr(7^7"7V) • 
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One then obtains 



tr(7''7"7'57^) 



tr(7''7"7''7^) 



where we have used 





(D.9) 



Pl{k - ViV + 2pi ■ (A; - pi){k - p,y -{k- pifpj = k\k- piy -{k- pifk^ , (D.IO) 

with k being one of the loop momenta. 

The contribution involving the one-loop vacuum polarization reads 



^ (47r)2-(4AA - 





-{P1+P2) 



\ 

(D.ll) 

Considering a factor —1 from the cancellation of the propagator connecting the chiral 
vertex to the two-loop self energy, we obtain from the D-algebra manipulations 

/ \- / \ 

1 




(47r) 



2XXp{ 



\ 




-m\-{y + \'))pi 




J 



(D.12) 

Let us conclude by mentioning a useful property that was used in appendix [El In 
the combination 

^ \ / (D.13) 
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the infrared divergence from the integrals involving the first leg is cancelled out. 
There are two other diagrams with non-trivial D-algebra and colour structure 




(D.14) 



Interestingly, these can be seen to be proportional to the very same integrals which 
appear in components [U]. The two diagrams are zero due to the vanishing of the 
one-loop triangle subdiagrams. 



E Cancellation of IR divergences 

In order to check the cancellation of the IR divergences, together with the contributions 
having both UV and IR divergences given in section HJ we have to include diagrams that 
have pure IR poles and would have been excluded by the UV finiteness conditions of 
subsection 14.11 The cancellation of IR divergences in the combination flD.13|) means that 



chiral and anti-chiral vertices with any number of legs and with external propagators 
are free of IR divergences from perturbative corrections. In the following we will hence 
attach propagators to the external fields of the diagrams that appear at four loops as 
quantum corrections of a chiral composite operator. This does not affect the UV poles, 
since the chiral self-energy is UV finite as demonstrated in appendix ID.2I 
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The contributions to the x(l) structure with only an IR divergence are given by 



Vr 



^ ^^MN{AMN - M2)/4irx(1) 

(4AA - A') -(-3 + 7 - ln47r)x(l 



Vr 



r45 



— > 



16 ■ ■ \ e 

-^^(4^^ - ^')(^4IR + /4UVIR - /4)X(1) 



AA 



(4AA-A2) -(-3 + 7-ln47r)x(i; 



Vr 



r46 



16 

(47r)^ MN 
AA / 2 

— (4AA - A^) -(-3 + 7 - ln47r)x(l) 
16 \ e 



{AMN - M^){hvYm -h + hhiK - K(/2)/2ir)x(1) 



Vr 



r37 



(47r)^ MA^ 
AA 



(4MiV - M2)(/4ir, + /2/21R - K(/2)/2ir)x(1) 



r38 



(4AA-A2) -(-3 + 7-ln47r)x(l) 
16 \ e 

{AtiY MN 



{AMN - M^){h - /4UVIR - /2/21R + K(/2)/2ir)x(1) 



AA / 9 

^(4AA-A2)(-(-3 + 7-ln47r)x(l) 
16 \ e 



(E.l) 

where we have given only the IR pole terms, and the UV sub divergences have been 
subtracted. 

We also have to consider the correction of the chiral propagator that is a neighbour 
of the fields interacting via x(l) 



V 



r3s 



i 



— > — 



{Any MN 



2 (8MiV-(M2 + iV2))(/2-K(/2))/2iRx(l) 



AA 



(8AA - (A^ + P)) ( Vs + 7 - ln47r)x(i; 



According to (]D.13p . one half of this contribution has to be taken into account, since 
the other half should cancel part of the IR divergence from an interaction of the isolated 
leg via a one- loop corrected gauge propagator with its neighbour to the right. Similar 
considerations hold also for the reflected diagram of V^^^i such that the total contribution 
of these diagrams to the IR divergence is |(1 + TVjVr^s- 

Further IR divergent contributions from self energy corrections of the three external 
and one internal line at the upper chiral vertex that forms cancel among respective 
diagrams in which two of these lines are interacting via one-loop corrected gauge prop- 
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agator. This is guaranteed by flD.lSP since in the considered propagators are attached 
to their external hnes. 

At this point a simple way to check the cancellation of the IR divergences is to sum 
up all the contribution containing them and check that the result is the same as if from 
the very beginning we had omitted all IR divergent diagrams, and had only considered 
V;uv and l^uv j^^ f^^^^ ^^e sum 

- (1 + 7^)(\/r35 + ^r44 + ^45 + V,A% + l^r37 + Kss) " 3^36 " ^(1 + ^)'^r3s 

= ^Miv( - GMNh^y,^ + ^{SMN - (M^ + N'))i3h + 2/42b6d))x(l) (E.3) 
turns out to be equal to 

-il+T^K'i^-W.'i^ , (E.4) 

which is the respective contribution of only the overall UV divergences from the diagrams 
with also an IR divergence to (14.101) . 

It is important to note that, besides the previously described check of the cancellation 
of the IR divergences, we have also performed the full computation of the range three 
contribution in the IR-safe r^-gauge described in [31]. The result turns out to be the 
same. 



F Double poles 

In this appendix we check explicitly the cancellation of the double poles in In Z. For that 
we need to consider diagrams which are responsible for interactions between magnons 
at odd and even sites which are proportional to chiral functions x(l; 2) and x(2, 3). We 
start by computing those contributions, and then we prove the complete cancellation of 
the double poles. 



F.l Odd- and even-site magnon interactions 

The relevant diagrams that couple the odd and even site magnons with each other are 
the following ones 



mixed 




16 



(-i + f)x(1.2) 



mixedl 




16 V4e2 26 



X(l,2) 



mixed2 




(47r)4 (MiV)2 (AA)2/ 1 U 

P 2 -'422btrylBc(iXll)2j — - + -jx(l,2j . 

(F.l) 
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In the sum of all contributions one has to consider the reflected diagrams. The second 
contribution acquires an additional factor of two due to two distinct positions for the 
vector vertices which are not mapped to each other under reflection. The result for the 
mixed renormalization constant read^ 

(AA)^ 1 

-^4, mixed (1 + 

(F.2) 

As expected [55], the 1/e pole is cancelled out such that at four loops there is no 
contribution to the dilatation operator that couples the magnons at odd and even sites. 



F.2 Double pole cancellation 

Summing up the contributions to the poles of the odd-site sector to the four- loop 
Z from (113!), (SSD, flCTD and (O}, we obtain 



,odd + ^4 ,mixed ) 1 



rS.odd 



4, mixed 



r4,odd 



+ Zr3,odd) I 4v 



[ \ 3) + x(3, 1)) + 2) - 



(F.3) 



In the definition of the dilatation operator, the logarithm guarantees that all higher 
order poles in e cancel out, such that \a.Z only contains simple ^ poles. Inserting f l3.2p . 
the expansion reads 



In Z = A^Za + A^ Z4 - -Z| + C(A 



(F.4) 



Let us now check the double pole cancellations in the A^ term. The two-loop contribution 
to the renormalization constant for operators of length L can be written as 



2L 



■^^(X(l) + X(2)) 



(F.5) 



where we have indicated the sum over the sites explicitly. It has an obvious decomposi- 
tion into two parts acting exclusively on even and on odd sites, respectively. The square 
of the above result can be decomposed as follows 



1 



22,dc 



22,5 



(F.6) 



^There is another contribution with identical prefactor that involves the chiral function x(2,3) that 
we associate to the even site sector. 
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The individual terms are given by 



22,dc 



22,5 



2L 

E 

j>i+3 
1 

-y 

i=l 



v 



+ 



+ 



+ 



+ 



(F.7) 



^ ^^M^iV^ Kihnxil, 3) + X(3, 1) + 2x(l, 2) - 2x(l)) 



2 fc4 
(AA)2 1 



(X(l,3) + x(3,l) + 2x(l,2)-2x(l)) 



16 2^2 

where the arrow denotes that in the final result we have considered the chiral functions 
with odd indices only and 2) and neglected the ones with only even indices and 
X(2,3). 

According to (IF.6p . the square of the two- loop contribution expands as 



(AA)^ 1 
16 2e' 



(X(l,3) + x(3,l) + 2x(l,2)-2x(l)) + ... 



(F.8) 



where we have neglected the chiral functions with only even arguments and x(2, 3). We 
have also disregarded the terms -Z22,dc which trivially cancel against four-loop diagrams 
that only contain double poles and hence become disconnected when the composite 
operator is removed. We have omitted to present these diagrams in the paper. 
Comparing equations (IF. 311 and (IF. 81) we finally find our desired result 



2 







(F.9) 



where we have considered that the discussion is identical for the neglected contributions 
with chiral functions with even arguments and x(2, 3). 
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